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Abstract 

We study global minimizers of a continuum Landau-De Gennes energy functional for nematic liquid 
crystals, in three-dimensional domains, subject to uniaxial boundary conditions. We analyze the physically 
relevant limit of small elastic constant and show that global minimizers converge strongly, in W^'^, to a 
global minimizer predicted by the Oseen-Frank theory for uniaxial nematic liquid crystals with constant 
order parameter. Moreover, the convergence is uniform in the interior of the domain, away from the 
singularities of the limiting Oseen-Frank global minimizer. We obtain results on the rate of convergence 
of the eigenvalues and the regularity of the eigenvectors of the Landau-De Gennes global minimizer. 

We also study the interplay between biaxiality and uniaxiality in Landau-De Gennes global energy 
minimizers and obtain estimates for various related quantities such as the biaxiality parameter and the 
size of admissible strongly biaxial regions. 

1 Introduction 

Nematic liquid crystals are an intermediate phase of matter between the commonly observed solid and liquid 
states of matter [13] . The constituent nematic molecules translate freely as in a conventional hquid but whilst 
flowing, tend to ahgn along certain locally preferred directions i.e. exhibit a certain degree of long-range 
orientational order. Nematic liquid crystals break the rotational symmetry of isotropic liquids; the resulting 
anisotropic properties make hquid crystals suitable for a wide range of physical applications and the subject 
of very interesting mathematical modelhng |18j . 

There are three main continuum theories for nematic liquid crystals [18j . The simplest mathematical 
theory for nematic liquid crystals is the Oseen-Frank theory [llj . The Oseen-Frank theory is restricted to 
uniaxial nematic liquid crystal materials (liquid crystal materials with a single preferred direction of molecular 
alignment) with constant degree of orientational order. The state of a uniaxial nematic liquid crystal is 
described by a unit-vector field, n(x) £ 5^, which represents the preferred direction of molecular alignment. 
In the simplest setting, the liquid crystal energy reduces to: 

^ofM = / n^,k(.x)ni,u{x) dx, (1) 
Jn 

the standard Dirichlet energy for vector-valued maps into the unit sphere. The equilibrium configurations 
(the physically observable configurations) correspond to minimizers of the J^iF-energy, subject to the imposed 
boundary conditions. In particular, the minimizers of J-^of are examples of S'^-valued harmonic maps |18|i32j. 
The Oseen-Frank theory has been extensively studied in the literature, see the review [2], and there are 
rigorous results on the existence, regularity and singularities of Oseen-Frank minimizers. 

The Oseen-Frank theory is limited in the sense that it can only account for point defects in liquid crystal 
systems but not the more complicated line and surface defects that are observed experimentally. A second, 
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more comprehensive theory is the continuum Ericksen theory [8]. The Ericksen theory is also restricted to 
uniaxial liquid crystal materials but can account for spatially varying orientational order i.e. the state of 
the liquid crystal is described by a pair. (s,n) e M x S*^. where s S M is a real scalar order parameter that 
measures the degree of orientational ordering and n represents the direction of preferred molecular alignment. 
In the simplest setting, the corresponding energy functional is given by 

jr^[s,n]= / s{x)^\\/n{x)\^ + k\\/s{x)\^ +Wo{s) dx (2) 
Jn 

where fc is a material-dependent elastic constant and Wq{s) is a bulk potential. The Ericksen theory is based 
on the premise that s vanishes wherever n has a singularity and this theory can account for all physically 
observable defects. 

However Ericksen recognizes that his theory is but a simplified description of a possibly more complex 
situation (see [8]): 

"There is the third possibility, that the three eigenvalues of Q are all distinct, giving what are called 
biaxial nematic configurations. Theories fitting MACMILLAN'S [11] format permit any of the three types of 
configurations to occur. Certainly it is not unreasonable to think that flows or other influences could convert 
a rather stable nematic configuration to one of the biaxial type, etc. I [19] am one of those who have argued 
that, near isotropic-nematic phase transitions, it should be quite easy to induce such changes. Accounting 
for such possibilities does add significant complications to the equations and the problems of analyzing them. 
Experimental information concerning the biaxial configurations is still quite slim and, for me, it is too early 
to think seriously about them. So, I will develop a theory representing a kind of compromise. " 

The most general continuum theory for nematic liquid crystals is the Landau-De Gennes theory [131 125j 
which can account for uniaxial and biaxial phases (biaxiality implies the existence of more than one preferred 
direction of molecular alignment). Indeed, this theory was one of the major reasons for awarding P.G. De 
Gennes a Nobel prize for physics in 1991. In the Landau-De Gennes framework, the state of a nematic liquid 
crystal is modelled by a symmetric, traceless 3x3 matrix Q £ M^^^, known as the Q-tensor order parameter. 
A nematic liquid crystal is said to be (a) isotropic when Q = 0, (b) uniaxial when the Q-tensor has two equal 
non-zero eigenvalues; a uniaxial Q-tcnsor can be written in the special form 

Q^s(^n(g)n-^Id^; seR\{Q}, nGS^ (3) 

and (c) biaxial when Q has three distinct eigenvalues; a biaxial Q-tensor can always be represented as follows 
(see Proposition [T]) 

Q = s (n(ii)n- ^Idj + r (m ®m - j s, r £ M; n, m G S'^ . (4) 



The Landau-De Gennes energy functional. J-lg[Q]: is a nonlinear integral functional of Q and its spatial 
derivatives. We work with the simplest form of J-lg[Q]i with Dirichlet boundary conditions, Qi, (refer to 
(fTTj) '). on three-dimensional domains O C R^. We take J-lg[Q] to be [27] 

:Flg[Q]= J^^\VQ\^{x) + fB{Q{x))dx (5) 

where fsiQ) is the bulk energy density that accounts for bulk effects, | VQp = j Qij,kQij,k is the elastic 
energy density that penalizes spatial inhomogeneities and L > is a material-dependent elastic constant. We 
take fsiQ) to be a quartic polynomial in the Q-tensor components, since this is the simplest form of /b(Q) 
that allows for multiple local minima and a first-order nematic-isotropic phase transition [131 132] . This form 
of fs (Q) has been widely- used in the literature and is defined as follows 

/b(Q) = ^i^^^tr (Q2) - ^tr (Q3) + £ (trQ^)^ 
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where a,b,c e R are material-dependent positive constants, T is the absolute temperature and T* is a 
characteristic liquid crystal temperature. We work in the low-temperature regime x < T* for which a{T — 
T*) < 0. Keeping this in mind, we recast the bulk energy density as follows: 

/b(Q) = -ytr (0^) _ ^tr (Q^) + ^ (trQ^f , (6) 

where a^6^,c^ € R+ are material-dependent and temperature-dependent positive constants. The equilibrium 
configurations (the physically observable configurations) then correspond to minimizcrs of Tlg[Q]^ subject to 
the imposed boundary conditions. 

In the first part of the paper, we study the the limit of vanishing elastic constant L —t Q for global 
mmmuzcrs, Q(^), of Tlg[Q]- This study is in the spirit of the asymptotics for minimizcrs of Ginzburg- 
Landau functionals for superconductors [3]. The limit L — > is a physically relevant limit since the elastic 
constant L is typically very small, of the order of 10^^^ Joule/metre. |27j . 

We define a limiting harmonic map Q(o) as follows 

where s+ is defined in (fTO|) . 71*^°) is a minimizer of the Oseen- Frank energy, Jx>fM in (III), subject to the fixed 
boundary condition n = rtf, G C°°{dQ., S'^) and Qb and are related as in pT|) . Our main results are: 

• There exists a sequence of global minimizcrs {Q^^>'^ } such that g(^^-) Q(°) strongly in the Sobolev 
space W^-'^{n,M.^). 

• The sequence {Q'-^''''} as above converges uniformly to as Lk 0, in the interior of fi, away from 
the (possible) sing ularities of Q("). 

• The bulk energy density, /s (Q'^*"^), converges uniformly to its minimum value away from the (possible) 
singularities of Q'-''-'; the uniform convergence of the bulk energy density holds in the interior and up to the 
boundary, away from the (possible) singularities of Q^^\ 

These results show that the predictions of the Oseen- Frank theory (described by the limiting map Q*-*^-*) 
and the Landau-De Gcnnes theory agree away from the singularities of Q^^\ The global minimizcrs, Q^^', 
are real analytic (see Proposition II 3p and have no singularities as such. However, one of the most intriguing 
features of nematic liquid crystals arc the optical 'defects' that appear in the Schlieren textures [13j . From 
a physical point of view, these defects are regions of rapid changes in the configurational properties of a 
nematic liquid crystal [13]. We conjecture that certain types of optical defects in Q^^*"-* (for small Lk), when 
they exist, may be localized near the analytic singularities of the limiting map Q*-"-*, since Q^-^''^ can have 
strong variations only near the singularities of Q''*''' (more precisely, the gradient, VQ^^''\ cannot be bounded 
independently of Lk on any set containing a singularity of Q*-"-*). There is existing literature on the location 
of singularities in harmonic maps [T] and this may allow one to predict the location of (optical) defects in a 
global Landau-De Gennes minimizer. 

Our convergence results analyze the limit of vanishing elastic constant L — > 0. Physical situations are 
modelled by small but non-zero values of the elastic constant L. Thus our convergence results show that for L 
sufficiently small, the limiting harmonic map Q^"' provides but a 'rough' description of Q'^^' i.e. Q^^' can be 
thought of as having a 'leading' uniaxial part plus a small biaxial perturbation, away from the singularities of 
(5*-°^. This small biaxial perturbation is of order 0(\/L) where L << 1 (see Section [5] for details). However, 
numerical simulations show that biaxiality may become prominent in the vicinity of defects [231 128] . In the 
second part of our paper, we study biaxiality and their role in global minimizcrs noting that biaxiality 

(if it exists) is one of the main differences between Q^^-* and the limiting approximation Q'^^\ More precisely, 
in Propositions fTTl and fT2l we obtain estimates for the size of the regions where Q^^*) can deviate significantly 
from Q'--^^ and on the size of admissible strongly biaxial regions in Q'^', in terms of the biaxiality parameter 
/3 (defined in (1^^ ) and the material-dependent constants. While Proposition [TT] mav be relevant to the 
properties of Q^^^ near the singular set of Q^*^', Proposition II 21 is relevant to the equilibrium properties away 
from the singular set of Q'^^K 
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Using a simple nearest- neighbour projection argument (see Corollary[T]) , wc show that the 'leading eigendi- 
rection', corresponding to the leading uniaxial part (see Section [5] for definitions) is smooth on any compact 
set K not containing any singularity of Q^'^\ Further, in Proposition [13 we also show that Q'^-' is either 
(a) uniaxial everywhere (except for possibly a set of measure zero where Q can be isotropic) or (b) Q*-^' is 
biaxial everywhere and can be uniaxial or isotropic only on sets of measure zero. It is known that as long 
as the number of distinct eigenvalues does not change, the eigenvectors of Q'^-' enjoy the same degree of 
regularity as Q^^' itself [26]. In Corollary [21 we show that the eigenvectors are necessarily smooth everywhere 
except for possibly a zero-measure set where the number of distinct eigenvalues changes and therefore, if the 
eigenvectors of Q^^-* suffer any discontinuities, these discontinuities must be localized on the uniaxial-biaxial, 
uniaxial-isotropic or biaxial-isotropic interfaces. This result may be relevant to the interpretation of optical 
data from experiments and wc hope to explore this connection in future work. 

Finally, we note that the Landau-De Gennes theory for uniaxial liquid crystal materials has strong analogies 
with the S-D version of the Ginzburg-Landau theory for superconductors [3] . The Ginzburg-Landau energy 
functional for a three-dimensional vector field, u : f2 ^ R'^, is typically of the form 

^gl[u]^ [ hyu\' + ^{l-\u\^)' dx (7) 

where e > is a very small parameter. The functional Tgl has been rigorously studied in the limit e — > 
which is analogous to the limit L ^ in our problem. The new mathematical complexities in the Landau-De 
Gennes theory for nematic liquid crystals come from the high dimensionality of the target space and also from 
the possibility of biaxiality in global energy minimizers. Future challenges include a better understanding of 
the qualitative properties of global minimizers for small but non-vanishing values of L, a better description 
of Q(^) near the singularities of the limiting harmonic map Q^^\ the regularity of the eigenvectors and 
eigenvalues, along with a deeper understanding of the appearance and role of biaxiality in global minimizers. 

The paper is organized as follows - in Section [51 we introduce the conventions and notations that are used 
in the rest of the paper. In Section [3 we state two representation formulae for Q-tensors that are useful for 
subsequent computations in later sections. In Section [^ we study the properties of global energy minimizers 
in the limit L ^ Q and prove the convergence results. In Section [3 we discuss the consequences of our 
convergence results and their relevance to the bulk energy density, the biaxiality parameter, the eigenvalues 
and the eigenvectors of a global Landau-De Gennes minimizer. In Section [SI we derive estimates for the bulk 
energy density, obtain bounds for the size of admissible strongly biaxial regions and discuss the interplay 
between biaxiality and uniaxiality in a global energy minimizer. 

2 Preliminaries 

We take our domain, Vi C M^, to be bounded and simply-connected with smooth boundary, d^. Let 5*0 C M'^'^^ 
denote the space of Q-tensors, i.e. 

So {Q e g,, ^ Q,, = 0} 

where we have used the Einstein summation convention; the Einstein convention will be assumed in the rest 
of the paper. The corresponding matrix norm is defined to be 

\Q\ v/t^Q^ - ^/g~Q~. 

As stated in the introduction, we take the bulk energy density term to be 

fB{Q) = -ytr (Q2) - ^tr (Q^) + ^ (tr(Q2))^ 

where a^,b'^,c^ G M arc material-dependent and temperature-dependent positive constants. One can readily 
verify that fsiQ) is bounded from below (see Proposition [HI [H]), and we define a non- negative bulk energy 
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density, fsiQ), that differs from fsiQ) by an additive constant as follows: 

f b{Q) ^ fsiQ) - mill fn iQ) . (8) 

It is clear that fsiQ) > for all Q E Sq and the set of minimizers of fsiQ) coincides with the set of minimizers 
for fsiQ)- In Proposition [HI we show that the function fsiQ) attains its minimum on the set of uniaxial 
Q-tensors with constant order parameter s+ as shown below 



= ^ Q e Qmin where 

n — 

3 



e 5o,Q = s+ ( n - i/d) ,n e §n (9) 



with 



= ^2 • (10) 



We work with Dirichlct boundary conditions, referred to as strong anchoring in the liquid crystal literature 
|13j . The boundary condition Qf, £ Qmin is smooth and is given by 

Qb = s+ (nt ®nb- hd^ , ub e {dn- 5^) . (11) 

We define our admissible space to be 

Aq = {Qe W^^'^ (r!; S-o) ; Q = Qh on 917, with as in (HH)} , (12) 
where W^-"^ (17; 5o) is the Sobolev space of square-integrablc Q-tcnsors with square-intcgrablc first derivatives 

1 /2 

I^. The corresponding W^^^^-norm is given by ||Q|]h'1.2(o) = (/n IQP + I^QP ■ In addition to the 

VF^'^-norm, we also use the L°°-norm in this paper, defined to be ||Q||l°=(ji) = ess sup2,gf2|Q(a;)| . 

We study global minimizers of a modified Landau-De Gennes energy functional, Fi^q\Q\^ in the admissible 
space Aq. The functional Flq\Q\ differs from Tlg\Q\ in © by an additive constant and is defined to be 



Flg{Q\^ I ^Q^,.k{^)Q^,.k{^) + fB{Q{x))dx. (13) 



For a fixed i > 0, let Q*-^^ denote a global minimizer of Flg[Q] in the admissible class, Aq. The existence of 
Q'^^ is immediate from the direct methods in the calculus of variations f^. The bulk energy density, fsiQ), 
is bounded from below, the energy density is convex in VQ and therefore, -Flg[Q] is weakly sequentially lower 
semi-continuous. Moreover, it is clear that Flg[Q] and J-lg[Q] have the same set of global minimizers for a 
fixed set of material-dependent and temperature-dependent constants ja'^, 6^, c^, L}. 

The global minimizer Q'^^ is a weak solution of the corresponding Euler-Lagrange equations |21j 

LAQi, = -a^Qy - (^Q.kQkj - ytr(g2)^ + c^Q.^tiiQ^) i,j = 1, 2, 3. (14) 

where the term 6^^tr((5^) is a Lagrange multiplier that enforces the tracelessness constraint. It follows from 
standard arguments in elliptic regularity that Q^"^' is actually a classical solution of and Q^^) is smooth 
and real analytic (see also Section 1^?^ . 

Finally, we introduce a "limiting uniaxial harmonic map" Q^^^ : fl Qmin', Q*"*^' is defined to be a global 
minimizer (not necessarily unique) of Flg[Q] in the restricted class, Aq D {Q : —^ Sq, Q{x) G Qmin s^.e. x € 
il}. Then Q^^-* is necessarily of the form 



Q 



(0) 



® - i/d^ , (15) 
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where is a global minimizcr of ToF[n] (see [2], [1]), 



|Vn°(x)|^dx= min / |Vn(x-)|^da; (16) 



in the admissible class An = G M^^'^ (^i; 5*^) ; n = rib on and rif, and are related as in (fTT|) . This 
"limiting harmonic" map Q*-"^ is therefore obtained from an energy minimizer, n", (not necessarily unique) 
within the Osccn-Frank theory for uniaxial nematic liquid crystals with constant order parameter (for more 
results about the relation between n'"^ and Q^"^ see [2])- It follows from standard results in harmonic maps 
[32] that Q'-*'-' has at most a finite number of isolated point singularities (points where has singularities). 
In the following sections wc will elaborate on the relation between Q^^-* and Q'^^\ 



3 Representation formulae for Q-tensors 

We have: 

Proposition 1 A matrix Q € Sq can be represented in the form 

Q ^ s{n®n-]-Id)+r{m(S)m-\-Id) (17) 
o o 

with n and m unit-length eigenvectors of Q, n ■ m — and 

< r < ^ or ^ < r < (18) 

The scalar order parameters r and s are piecewise linear combinations of the eigenvalues of Q. 

Proof. From the spectral decomposition theorem we have 

Q = AiTii <g)ni + X2n2 ® n2 + Xsn^ (g) (19) 

where Ai, A2, A3 are eigenvalues of Q and rti, rt2, are the corresponding unit eigenvectors, pairwisc perpen- 
dicular. We have / = X]?=i ^ ^^'^ ^^^^ tracelessness condition imphes that Ai + A2 + A3 = 0. Thus 

Q = Aini (g) ni + A2n2 ® n2 — (Ai + A2)(/ — ni (g) ni — 712 ® ^2) 

We consider six regions Rf ,i = 1, . . . , 6 in the (Ai, A2) - plane which cover exactly half of the whole plane. 
This corresponds to the representation ([T7)) with < r < | . The other half of the plane is covered by the 
regions R~ , i = 1, . . . , 6 , (which arc obtained by reflecting through the origin (0, 0)) and the regions R~ 
correspond to the representation (|17p .with r, s < 0. 

We let R^ = {(Ai,A2) G -2Ai < A2, Ai < 0}. In this case r 2Ai + A2 and s 2A2 + Ai with 

def def 

n = 712, m ~ ni. One can directly verify that for r, s thus defined, we have 

r = 2Ai+A2<^ = A2 + y. 

Interchanging Ai with A2 in the definition of r and s and m with n, we obtain the region R.^ = {(Ai, A2); A2 > 
-Ai/2;A2<0}. 

Let R^ = {(Ai, A2) e M^^ A2 < 0, A2 > Ai}. Taking r "^^ A2 - Ai, s -2Ai - A2, n =^ 713, tti n2, one 
can check that 

r = A2 - Ai < - = -Ai - y. 
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The region is obtained from interchanging Ai and A2. 

We have R+ = {(Ai, A2) e M.'^, Xi < 0, -2Ai > A2 > -Ai} with r =^ -2Ai - A2, s =^ A2 - Ai, 71 = 712 and 

def 

m = n^. Again, it is straightforward to check that 

Interchanging Ai with A2, we obtain the region Rg . 

FinaUy the remaining half of the ( Ai , A2 )-plane is covered by the regions (obtained from R^ by changing 
the signs of the inequaUties and keeping the definitions of r and s unchanged). For example, is defined to 
be 

R^ = {(Ai,A2) e m2;Ai > 0, 2Ai < -A2} 
with r = 2Ai + A2 and s = 2X2 + Ai. One can then directly check that 

s 

- < r < 0. 
2 ~ ~ 

The remaining five regions for i = 2 . . .6 can be defined analogously. □ 

Remark 1 The representation formula JiTp is known in the literature \25^ . In Proposition^^ we state that 
it suffices to consider the two cases given by il8\) : we have not found references for this fact. 

In Proposition [2 we state a second representation formula for admissible Q £ Sq and its relation to the 
representation formula (|17p . The representation formula (|20|1 is known in the literature [22j and will be used 
in Section m For reader's convenience we provide a quick proof. 

Proposition 2 (A second representation formula) A matrix Q Cz Sq can be represented as: 

Q^s(^®n~^Id^+R{m®m-p®p) (20) 

The vectors n, m and p are unit-length and pairwise perpendicular eigenvectors of Q with corresponding 
eigenvalues Ai, A2,A3. The scalar order parameters S and R are given by 

S = 3^ i? = i (2A2 + Ai) . (21) 

Proof. We have the spectral decomposition of Q, namely 

Q = Ain ® n + A2TO (g) m + X^p (g) p 
with n,m,p pairwise perpendicular unit-length eigenvectors of Q and 

Id = n'S)n + m(g)m + piSip. 

Combining the last two relations and taking S ~ 3^, R~\ (2A2 + Ai) we obtain the claim. □ 

4 The limiting harmonic map 

4.1 The uniform convergence in the interior 

Firstly, we recall that for a Q e Sq the biaxiality parameter (3{Q) (see for instance [23]) is defined to be 

/3(Q) - 1 - 6-^^^ (22) 
(trQ2)3 



The significance of P{Q) as a measure of biaxiality is due to the following 
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Lemma 1 (i) The biaxiality parameter f3{Q) e [0, 1] and (3{Q) ~ if and only if Q is purely uniaxial i.e. if 
Q is of the form, Q = s {n® n — \ld^ for some s G M, n G S^. (ii) The biaxiality parameter, P{Q), can be 
bounded in terms of the ratio ^, where (s, r) are the scalar order parameters in Proposition^^ . These hounds 
are given by 

^ <-<ifl + \/l-V^V (23) 



s - 2 



Equivalently, 



where (S, R) are the order parameters in Proposition\^ Further (3[Q) = 1 if and only if r ^ ^ or if and only 
if ^ = i. (iii)For an arbitrary Q G 5*0, we have that 



^/6V "^J ~ \ 2 

Proof: The proof of Lemma [T] is deferred to the Appendix. □ 
The next proposition gives us apriori L°° bounds, independent of L. 

Proposition 3 Let <ZK^ be a bounded and simply- connected open set with smooth boundary. Let Q*-^-* be 
a global minimizer of the Landau-De Gennes energy functional \1SI\) , in the space 111 



LAQ,j = -a^Q,j - b^ Q,pQpj - -trQ'^S.j ] + (trQ^) Q^^, (27) 



Then 

||Q^^)||L^(n) < (26) 

where S-|_ is defined in ilU\} . 

Proof. Proposition [3] has been proven in |21] : we reproduce the proof here for completeness. 
The proof proceeds by contradiction. In the foUowing we drop the superscript L for convenience. We 
assume that there exists a point x* G where |Q| attains its maximum and |Q(x*)| > ^J%^s+. On dVt, 

\Q\ — \J\sj^ by our choice of the boundary condition (note that if Q G Qmin then |Q| = y^.s+). If Q is a 

global minimizer of F[Q] then Q is a classical solution (see Section 16.21 for regularity) of the Euler-Lagrange 
equations 

1 

3^ 

Since the function jQp : f] ^ M must attain its maximum at x* G f2, we necessarily have that 

AQ|g|^)(x*)<0 (28) 
We multiply both sides of ([7f|) by Qij and obtain 

L A QlQP^ = -a^trQ^ _ b\^Q^ + (trQ^)^ + L|VQ|^ (29) 

We note that 

- aMQ^) - bMQ') + (tr(Q2))' > /(|Q|) (30) 

where 

f{\Q\) = -a^\Q\'-^\Q\' + c'\Q\\ (31) 
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since tr(Q'^) < -l^L from ([25)1 . One can readily verify that 



/(|Q|)>0 for|Q|>y^s+ (32) 
which together with ([29]) and (|30|) imply that 

AQ|Qp)(a;)>0 (33) 

for all interior points x G il, where |Q(x)| > y^.s+. This contradicts ([25)) and thus gives the conclusion. □ 

In what follows, let eL{Q{x)) denote the energy density eL{Q{x)) ^|VQp + We consider the 

normalized energy on balls B{x^ f) d ^ ^ {y G ^;\x — y\ < r} 

J'{Q,x,r)''=^ - [ eL{Q{x))dx = - [ + ^ |VQ|' dx. (34) 

JB(x,r) JB[x,r) ^ ^ 

We have: 

Lemma 2 (Monotonicity lemma) Let Q'^^ 6e a global minimizer of Flq[Q] in S13\) . Then 

J'{Q^^\x,r) < J'{Q^^\x,R),yx £n,r < R, so that B{x, R) C fl (35) 



Proof. The proof follows a standard pattern (see for instance and is a consequence of the Pohozaev 
identity. We assume, without loss of generality, that x = and < i? < d{0,dft), where d denotes the 
Euclidean distance. Since Q^^^ is a global energy minimizer, it is a classical solution (see Section 16.21 for 
regularity) of the system ([Tl)) : 

In pop and in what follows, we drop the superscript L for convenience. 

We multiply ([5S)) by Xk-dkQij, sum over repeated indices and integrate over B{0, R) to obtain the following 



dfsiQ) ,2% 
9Qy 3 



tr(Q2) 



(36) 



= / Q,,,o(.t) • Xk ■ d,Q,,ix) - l^M^M . . d,Q,j{x) - h^%r{Q\x)) ■ Xk ■ dkQrjix) dx 
Jb{o,r) ^ oQij ^ o 



1 dfBiQ jx)) 

B{aM) Jb{o,B) ^ "'^ 



Qij,u{x) ■ Xk ■ dkQijix) dx - I -j- — -Xk-dkQijdx (37) 



/ // 
where we have used the tracelessness condition Qu ~ 0. 
Integrating by parts, we have that: 

/ = / Qij,ii{x)xkdkQij{x)dx 

JB(0,R} 



i].i{5ikQij,k{x) + XkQij,ki{x))dx + I Qi],iXkQij,k-^ dx 



xi_ 

B(0,R) JdB(0,R) ^ 

1 



ij,i{x)Qij^i{x) dx + 3 / T:Qij,i{^)Qtj,i{x)dx 



B(O..R) Jb{om) 2 

,,(x)Q,,,(^) Xk ■ Xk r {Q.M-Xkr (38) 



dB{0,R) 2 R JaB{0,R) R 
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JB{0,R) ^ OLjij L Jb{O.R) 

^ ^ MQix)) dx + y f fsiQix)) ■ dx (39) 



L ■'b{0,R) L JdB{O.R) ^ 



Hence (l37l) becomes: 



}ij,iix)Qij,i{x) fsiQix)) f Qij,k{x)Qij,k{x) , fsiQix)) 



B{Q,R) ^ ^ JdB{Q,R) 2 L 

I {Q,,,k{x)-Xkfdx + 2 I MMldx (40) 

■ldB(Q.R} JB{Q,R) ^ 



We have 



Q.,,i{x)Q^3,i{^) I fB{Q{x)) A ^ 1 f Qij,l{x)-Q,,,l{x) ^ fB{Q{x)) 



dR \ R .1 B(o,R) 2 L J R^ Jb(q.r} 



R 

The right-hand side of (|^T|) is positive from (|in)) and hence the conclusion. □ 



1 f Qi].i{x) ■ QijA^) _^ fBiQjx)) ^^^^ 

dB{0,R) 2 L 



Lemma 3 ( W^''^ convergence to harmonic maps) Let f2 C M'^ he a simply- connected bounded open set with 
smooth boundary. Let Q'^' be a global minimizer of F^gIQ] *^ the admissible class Aq defined in il2]) . Then 
there exists a sequence Lk ^ so that Q^^*-' Q^") strongly in W^'^ {Q; Sq) , where Q'--^^ is the limiting 
harmonic map defined in il5\) . 

Proof. Our proof foUows closely, up to a point, the ideas of Proposition 1 in [S]. Firstly, we note that the 
limiting harmonic map Q^^^ belongs to our admissible space Aq and since Q^^^{x) € Qmin, a-G- x £ fl (see 
Section ^ we have that fs {Q^°Hx)) = a.e. x & n. Therefore 



/ dx< [ IqI^I{x)QI^I{x) + }-UQ^^){x))dx < [ iQS(x)Q(y,(x)dz (42) 

The Q'-^^'s are subject to the same boundary condition, Qb, for all L. Therefore shows that the 
M^^'^-norms of the (5*-^'''s are bounded uniformly in L. Hence there exists a weakly-convergent subsequence 
Q'-^''^ such that Q^^*") ^ Q^^^ in W^'^, for some Q^^^ g Aq as Lk 0. Using the lower semicontinuity of the 
p^i,2 ]2orm with respect to the weak convergence, we have that 

/ \VQ^^\x)\'' dx < [ \VQ^''\x)\^ dx (43) 
Jn Jn 

Relation g2]) shows that fBiQ^^-Hx)) dx < Lk JnQif,ki^)Qif,ki^^ ^^"^ ^^^^^ / hiQ^^'-Kx)) dx 
as Lfc 0. Taking into account that /b((5) > 0,V(3 G 5*0 we have that, on a subsequence Lkj, 
fB{Q^^''^\x)) — *■ for almost all x & fl. From Proposition [HI we know that fsiQ) = if and only if 
Q G Qmin i-c. if Q = s+ (n (g) n — \ld) for n e S^. On the other hand, the sequence Q^^*^) converges weakly 
in W'^''^ and, on a subsequence, strongly in L^ to Q^^'. Therefore, the weak limit Q'^^ is of the form 

Q^^\x) = s+ fn'-^\x) ® n^^\x) - ^Ldj , n^^\x) G S^ a.e.xe n (44) 
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It was proved in (see also [4]) that if Q''^'' & M^^'^ and the domain is simply-connected, we can 
assume, without loss of generality, that n^^^ S W^'^{^, S^) and its trace is Tif,. Then dBl) implies \\7 Q'-^\x)\'^ = 
2s^|Vn(^'(x)p for a.e. x e fl. Also, recalling the definition of Q'") from Section [2]we have \VQ^°\x)\'^ = 
24|Vn(°)(a;)|2 for a.e. x e fl. 

Combining ([^5]) with and the observations in the previous paragraph, we obtain J^^ |Vri(^^(a;)p dx = 
^ |V7i(0)(x)|2dx and |Vg(i)(a;)|2 = ^ |VQ(°)(x)|2 dx. Then: 

/ |VQ(°)(a;)pda; < liminf / IVg^^'^^x)! < limsup / |VQ^^'"^Hx)p dx < / |VQ(°)(x)P dx, 

which demonstrates that limLfc.-»o ||VQ'-^'"j''||l2 = ||VQ'°-'j|i2. This together with the weak convergence 
Q^^kj) _^ q{o) gufgces to show the strong convergence q'^'^j-' Q^"' in W^-"^ . □ 
The following has an elementary proof, that will be omitted: 

Lemma 4 The junction /b'-Sq^ R+ is locally Lipschitz. 

We can now prove the uniform convergence of the bulk energy density in the interior, away from the 
singularities of the limiting harmonic map Q^^\ 

Proposition 4 Let d be a simply-connected bounded open set with smooth boundary. Let Q^^-* € 
W^'^{fl, So) denote a global minimizer of Flq\Q] in the admissible class Aq. Assume that we have a sequence 
{Q^^-^keN so that Q^^k) ^ g(o) j„ w^''^{n, So) as Lk 0. 

For any compact K <Z il such that Q^*^' has no singularity in K we have 

Vnn fB{Q^^''\x))=Q x G A' (45) 

Lfc^O 

and the limit is uniform on K . 

Proof. Lemma |3] shows that the strong limit Q'^"-' is a limiting harmonic map, as defined in Section [51 
(5^°) = s+(n^'''(x) ® n*^°-'(x) — where n''''(x) € W^^'^(r2,S^) a global energy minimizer of the harmonic 
map problem, subject to the boundary condition n = nf, on 917. 

Let aj^k = /s((3'-^'°-'(xo)), for xq G K an arbitrary point. Proposition [3]and Lemma pimply that there 
exists a constant (3 (independent of xq) so that 

l/i3(Q'^'(^)) - hiQ'^'^Hv))] < - Q'^'(y)l (46) 

for any x,y ^ ft, L > 0. 
We then have 

< ./b(Q'^''^(x)) +/3|Q(^'=)(x) - g(^''^(xo)| 
< MQ^'-'-Hx)) + PWVQ^^-^h^^K'jloo ~ xol < fB{Q^'^'\x)) + -^|x - xo|,Vx e K' (47) 

where K' C is a. compact neighborhood of K to be precisely defined later. In the last relation above we 
use Lemma A.l from [3] and the apriori bound given by Proposition [3l For reader's convenience we recall 
that Lemma A.l in [3] states that if u is a scalar-valued function such that — Au = / on C M" then 

|Vm(x)P < C (^||/||L~(0)llu||L~(n) + dist'-'("l.3a) ll"Hi°°(o)) where C is a constant that depends on n only. In 
our case the constant C depends on the dimension, n = 3, on a^,b^,c^ and on the distance sup^^^ d(y, 3f2) 
only. 
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From (|47|) we have that 

a^. - < fB{Q^'"'\x)),yx e K\ \x - xo\ < Pk (48) 

We argue similarly as in [5] and divide by Lk and integrate over Bp^{xQ) to obtain: 

^(a..-%^)</ M^,, (49) 



Take an arbitrary e > 0. Recall that isT is a compact set that does not contain singularities of Q'"'- 
Then there exists a larger compact set K' , so that K C K' , that does not contain singularities either, and a 
constant Ck' such that \VQ^°\x)\'^ < Ck', Vx G K'. For Rq small enough, with Rq < dist{K,dn) and such 
that B{xo,Ro) c K' yxo G K we have 

If W^,,<f^C.^«.<|.V.,.^ (50) 

We fix an i?o as before. As Q^^'''> Q'"' in W^'^, we have that there exists an Zq > so that: 



^oJBn,,{xo) 2 RoJbr^(xo) 2 3 

The arguments in [3] fail to work in our case as we have a three dimensional domain, unlike in the quoted 
paper, where the domain is two dimensional. In our case, using the monotonicity formula from Lemma 2 and 
taking pk < Rq we obtain: 

/.(Q(-^)(x)) < ^ / + M^^,. < + (52) 



lB,^{xo) ^oJBn,,{xo) 2 Lk V 3 3, 

for Lk < Li with Li small enough so that -s- fr, , ^ -fatQ' '^H^)) < % (note that there exists such an Zi 
as the proof of Lemma 3 shows that J^^ ^'^'"'^l^ ^'^ ~ "^^^^ ^ ~^ 
We take pk = "''^^ ■ Then, from (gH) and ^ we obtain 

for L;; < min{Zo, Li}. As e > is arbitrary and the estimate on a^^ ~ /-b(Q'"^'°H^o))j xo & is obtained in 
a manner independent of xq, we have the claimed result. □ 

We also need the following 
Lemma 5 There exists Eq > so that: 



2 



i,j=l \ / 

VQ S So such that \Q — s+(n (E) n ~ — /d)| < eo, for some n £ (53) 

o 

where s+ = '''^^4^2^'^°^'^^ a'^c^ i^e constant C is independent of Q, but depends on a^,5^,c^. 
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Proof. Recall from Proposition [HI [21] that fsiQ) > and fsiQ) = ^ Q = s+{n ®n - \ld) with 

Let the eigenvalues of Q be a;, y, —a; — ?/. We define F{x, y) — a^(x^ + + a;?/) + h'^xy{x + j/) + c^(a;-^ + 
2/2 + xyf and D =^ min(,,^)eK2 F(a:, y). Then F(.t, y) - = /b(Q). 

Then F = only at three pairs (x, i/) namely (— ^^): i^^^ ^^"^ (^^' ~^)- 
On the other hand we have 



i,j = l \ / 

+c4(tr(Q2))3 + 2a262tr(g3) - 262c2tr(Q2)tr(Q3) (54) 



(where we used the identity tr(Q'*) = llli^I-L^ valid for a traceless symmetric 3x3 matrix) 



by 



2 

If we denote h{Q) = Y.h=i {^^^ + ^^^-^(Q^)^ we have h{Q) = H{x,y) where iJ : ^ M is given 



, r hi 

H(x, y) = 2a^(x^ + y^ + xy) + 4( 2a^c^){x'^ + y^ + xyf + Sc^fa-^ + y'^ + xyf 

6 

+ 12b^c^xy{x + y)(a;2 + y'^ + xy) — 6a^b^xy{x + y) 
We claim that there exist ei,e2,e3 > so that 

ipix^y) < H{x,y) < CF[x,y)y{x,y) e B,^ (-^, -^), B,, (-^, 2^), 5,3 (2^, (55) 

which gives the conclusion. 

We prove the inequality ([55]) only for {x,y) G Bs^{—^,~^); the other two cases can be dealt with 
similarly. 

Careful computations show: 

s+ s+ dH s+ s+ OH s+ s+ 



dxdy^ 3' 3^" ^ 24c4 



^ 3' 3^ dx^ 3' 3^ 3' 3^ 



a?;2 ^ 3 ' 3 ' ' 3 ' 3 ^ 4^2 



92^ , s+ s+ 



(-^, = + 12a'c2 + 52^54 + 24a2c2) 



2 



9x9?/ 3 ' 3 
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Let (a;o,yo) = We have 



F{x,y) Fi{x,y) + Rp.{x,y) 

where Hi{x,y) = {x - xo)'^^{xo,yo) + 2{x - xo){y - yo)-§^{xo,yo) + {y - ya)'^^{xo,yo) and Fi{x,y) = 
{x - xo)'^^{xo,yo) + 2{x - xo)(y - yo)^^{xo,yo) + {y - yo)^|p^(a;o, yo) with Rh, Rp the remainders in 
the Taylor expansions around (x'o,yo)- 

From the definition of Taylor expansions, we have that there exists eo > so that on B^^ {xq, yo) we have 
\RH{x,y)\ < ^Hi{x,y) and \Rp{x,y)\ < ^Fi{x,y)y{x,y) e Be^{xo,yo) (56) 

On the other hand we have 

^i(^'2^)^77;CTT^ ^ ^i(^'y) ^ Fi(x,y)8(54 + 6aV)V(x,2/) G (57) 
hence, combining ([56|) and ([57|) . we get: 

24(^,4+^^2^2) =^ ^(^' y) ^ ^(^' y)24(^' + Qa'c^^ix, y) £ S,, {-'-j-'-j) (58) 
which yields claim ((55)) for (x, y) G -Bsj ~^)- "^^^"-^ other two cases can be analyzed analogously. □ 

We continue by proving a Bochner-type inequality that is crucial for the derivation of uniform (in L) 
Lipschitz bounds, away from the singularities of the limiting harmonic map. This type of inequalities were 
first used (to the best of our knowledge) in the context of harmonic maps (see [29] and the references there) 
and later adapted to other, more complicated contexts (see for instance [5|). The main difficulty in the proof 
of Proposition [5] (to follow) is the derivation of the next lemma. 

Lemma 6 There exists Eq > and a constant C > 0, independent of L, so that for Q^^^ a global minimizer 
of Flg[Q] the admissible space Aq, we have 

-Aei(Q(^))(a;) <Cei(Q(^)(x)) (59) 

provided there exists a ball Bp(^x-j{x) for some p{x) > such that \Q^^\y) — s+ {jn{y) ® rn{y) — ^Id) \ < Eq 
with m{y) e S^, for all y G Bp(r^'^{x). 

Proof. In the following we drop the superscript L for convenience. We have: 



Qij.kQijjk 



-A{Q 



< -d. 



LdQ^ 



3L 



.tr(g2) 



Qij, 



ij.k)Qij,k * 

1 dfs 



< 



LdQ, 



'ij,k 



(60) 
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On the other hand: 




L dQii ^"^ ^ " ' 3L 




<-dkU§^iQix))]Q..,,k (61) 



We take ei > a small number, to be made precise later. For any such ei we can pick £o > small enough 
so that if the eigenvalues of Q{x) are (A, /i, —A — ^) then one of the three numbers (A + + (// + + 



iX + fi + 2^)^ (A + + {fi- 24±)2 + iX + fi- (A - 24±)2 + + £±)2 + (A + /X - is less' than 
or equal to ei (this can be done because the eigenvalues are continuous functions of matrices, [17], and the 
matrix s+{n(^ n — ^Id) has eigenvalues and 2-^). Note moreover that we need to choose Eq to be 

smaller than the choice (of Sq) in Lemma Oas we will need to use that lemma in the remainder of this proof. 

For the matrix Q{x), let us denote its eigenvectors by ni{x) , n2{x) , n3{x) and let Ai(a;), A2(a;) , A3(a;) — 
— Ai(x) — X2{x) denote the corresponding eigenvalues. From the prececding discussion, we can, without loss 
of generality, assume that 



(Ai + ^)^ + (A2 + ^)^ 



(Ai+A2 + 2^)2 <ei 



(62) 



We define the matrix 

„^ de_f s+ 



ni{x) (g) ni{x) — '^n2{x) ® 712(2;) H — ^n^lx) ® 713(2;) 



3 ^ ' ' ' 3 ' ' ' ' 3 
(Note that there exists a ■m{x) G §^ so that = s+(m(x) (g) m{x) ~ ^Id)). 

Taking into account and the fact that Q{x) and have the same eigenvectors, we have 



tr(Q(.T) 



r = (Ai + fr + (A2 + fr + (Ai + a2 + 2^^)^ < s. 



(63) 



Using the of Taylor expansion of ^ qq^. qq — {Q{x)) around wc obtain: 



^'^^ (Q(x)) = 1_^!^(Q-) + fJ''^^ (QHQpgix) - Q^J +7^^^"™(g^Q(a;)) (64) 



where TZ''^™'" {Q'^ , Q{x)) is the remainder. 
From we have: 
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^ij,k — J- ' 



1 d^fB 



~Qrfin,kQ 



1 d^fE 



L dQ.jdQ 



1 



L dQijdQrnndQ 



PI 



mn 



<0 



1 

d'fB 



^2 ^,l,ra,.=l ^ dQ,,dQ„,ndQ. 



7^'^'""(g(x),Q-)g,,■fcQ™„,fc< 



-(Q^) (gp,(.^)-Q^/ 



^ ^ (7^^^"")'(Q(x),Q-) + i|VQ|4< 



■i,j,m,n— 1 



2,j,m,n— 1 



< ^tr(Q(x)-g-)2 + i|VQ|4 



(65) 



where < 5 < 1 and Co, Co, Ci are independent of L and x. For the first term in the second line above we 
use the fact that the Hessian matrix of a function fsiQ) is non-negative definite at a global minimum (which 
holds true in our case as well, as one can easily check, even though we have fsiQ) restricted to the linear 
space So)- 

Let us recall (from the proof of the previous lemma) the definitions of F and F. Then, for a matrix Q Cz Sq 
with eigenvalues (Ai, A2, — Ai — A2) we have 



/b(Q) = F(Ai,A2) 

We claim that for ei > small enough there exists C2 independent of L, Ai, A2 so that 



(66) 



C2 ((Ai + ^ f + (A2 + ^)2 + (Ai + A2 + 2^)2) < F(Ai, A2) 



for all (Ai, A2) so that(Ai + -jY + (A2 + -fy + (Ai + A2 + 2^y < ei. (67) 



Careful computations show: 



^ 3 ' 3 ^ d\i^ 3 ' 3 ^ 9A2 ^ 3 ' 3 ' 



d^F , s+ s+, d^F , s 



9A2 ' 3 ' 3 ' dXl' 3 ' 3 ' 4c2 



-( — 5"' — T) = 3a 



dXidM ' 3 ' 3 
Using a Taylor expansion around (Ai, A2) = {—^, ~l3") have 
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^(Ai, A2) = ^ + 12a'c2 + 6V&^ + 24a2c2) [(Ai + '^f + (A2 + ^f] + 

+3a2(Ai + ^)(A2 + y) + i?(Ai, A2) > 
> ^ 1^ (^^^ + I2a'c' + b'Vb^TM^) [(Ai + + (A2 + + 3a'{X, + ^)(A2 + ^)| (68) 

where R{Xi , A2) is the remainder in the Taylor expansion, and the inequahty holds provided that the remainder 
i? is small enough. We choose ei > to be small enough so that if (Ai + ^)^ + (A2 + 3-)^ + (Ai+A2+24^)^ < ei 
then R is small enough and the inequality above holds. 

As the quadratic form {b^ + Ua^c^ + b^Vb^T2A^) [(Ai + ^ )2 + (A2 + + fa^ (Ai + ^) (A2 + ^) 
is positive definite, there exists a C2 > 0, depending only on a^,b^ and such that 

l[^{b' + 12aV + bWb^ + 24a^c^) [{X, + 'ff + {X, + 'ff] + 3a'{X, + '-±)(X, + ^)} > 

> C2 (^(Ai + + (A2 + + (Ai + A2 + V(Ai, A2) e 

Combining this last inequality with we obtain the claim ((57|) . 

The relation §7i) together with jMl) and §^ show that tr(Q(a;) - Q^)^ < CsfsiQix)) for some C3 
independent of L and x, which combined with (|65p shows 

with C4 a constant independent of L and x and any S > 0. This last inequality together with ([HO]) and ((^T|) 
show: 

Taking into account Lemma O and choosing iS small enough (depending only on C4 and the constant C 
from Lemma^ we can absorb the term ^^fsiQ) on the right hand side into the left hand side and obtain 

-Aei < i|VQ|^ 

giving the desired conclusion. □ 

Lemma 7 Let fl G M.^ be a simply-connected bounded open set with smooth boundary. Let Q^^^^ <E W^'^{fl, Sq) 
be a sequence of global minimizers for the energy Flg[Q] in the admissible space Aq. Assume that as Lk 
we have Q^^"^ Q^°^ in W^^'^{n,So). 

Let K G be a compact set which contains no singularity of Q^^\ There exists Ci > 0, C2 > 0, Zq > 0(^ 
all constants independent of Lk) so that if for a G K,0 < r < d{a, dK) we have 

- / eiJQ(^'=)(x))da;<Ci 

then 

r^ sup eL,(Q(^''))<C2. 

B.(a) 

for all Lk < Lq. 
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Proof. Taking into account our assumptions on the sequence (Q^^*"') ^.^pj, Proposition |4] shows that for 
any given £q smaller than Eq in LemmaE] and also smaller than the Eq in Lemma^ we have that there exists 
a Lq so that for < Lq we have 

||g(ifc)(^) - .s+ (^n{x) (g) n{x) - ^/d^ || < Eq, V.t € K, for some n{x) £ (69) 

Wc continue reasoning similarly as in |29j . We fix an arbitrary Lk < Lq and an a G O and take a r > so 
that < r < mm{d{a, dil) , d{a, K)} . Wc let ri > and Xi e -Bri(a) be such that 

max (|r-s)2max eL,(Q(^'')) = (rr-nf max eL,(Q(^'=)) = (|r - r^^e^, (Q^^'-^ (xi)) 

0<s<2r O B,(a) d -B^^ (a) 6 

Define e^^'^ ''^ maxB^^(a) et,(Q(^^)). Then: 



max eL,(Q(^'=))< max eL.iQ'^^"^) 

B -2/3 r-T^ (xi) T + (g) 

2 2 

(2/3-r-ri)2maxB (al ei^. (Q^^''-') ,r n n \ 

- (2/3T-(2/3.r + ri)/2)2 ~ ^ i^'fw ^^-^ ^ ~ ^''^ ^'"^^ 

where for the first inequality we use the fact that Sg/sr-rjj {xi) C i? (2/3 r+ri) (a) and for the second inequality, 

2 2 

we use the definition of ri . 



Let rs = (2/3-'- n)V^ ^^^^^ = Q^^") + j ■ Wc let Lu = e^^"^ Lu and then 



1 X 

Equation fTUl) then implies 



max e, = max ^^-^^f'^'^ < 4 

1 ^ 

where the equality above follows from the definition of r2 and R^^''^ and the inequality above follows from 
equation (|70p . Thus, wc have 

max e£ji?(^'=)) < 4, ez, (i?(^'=))(0) = 1 (71) 

^'■2(0) 

where R'^^'''^ satisfies the following system of elliptic PDEs 

LkR^:i = -a'R^^"^ - [R^t^R^^,"^ - ^tr((i?(^'^))2)j + c2i?,(f )tr((i?(^'=))2) (72) 
Wc now claim that 

r2 < 1 (73) 

It is clear that ^2 < 1 implies the conclusion. Let us assume for contradiction that r2 > 1. Then wc claim 
that there exists a constant C > , independent of L^, so that 
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1<C (74) 



The matrix i?'^-^'') satisfies the system ((7^ (which is the rescaled version of ((Ti)) ); using relation ((H^ and 
the definition of i?*^-^*") as well as the fact that r2 > 1, we can apply Lemma [5] to (i?'^-^'=') and obtain 

-Aei^iR^^-Hx)) <Cel^{R'^^''\x)) < 4Ce£ji?(^'=)(x)),Vx G Bi(0) 

Combining ([7T|) and the Harnack inequality (see for instance [HI], Ch.l4, Thm. 9.3) along with the above 
relation we obtain ([7i|) . 

We have 

/ ei^{R^ '>{x))dx<- + dx^ 



2/3 • r - ri 



e,.JQ(^'^)(x))dx<- / e^AQ^^^Hx)) 



< - / eL,iQ'^''''Hx))dx<3Ci (75) 



r 



iB^a) 

where for the first inequality we use the monotonicity inequality (Lemma [5]) and the assumption that r2 > 1 
(note that the equation satisfied by i?'^*"^ , equation ([7^ is the same as the equation satisfied by Q^^''\ up 
to a different elastic constant, hence the use of Lemma [2] here is justified). For the equality in relation ([75|) 
we use the change of variables y ^ xi -\ — ^== and use the relation: e£^(i?''^'=-') — (2^) eLfc(Q^'^ '"'')■ For the 



second inequality in (|75p we use the monotonicity inequality and the fact that < r. Por the third 

inequality in ([75)1 we use the fact that Bj,/^{xi) C Br{a) since l^i ~ a\ < ri < |r. The last step in ([75]) follows 
from the hypothesis of the Lemma. 

Choosing Ci small enough we reach a contradiction with (j74p which in turn implies that r2 < 1 and hence 
the conclusion. □ 

We can now prove the uniform convergence of Q^^'''^ away from singularities of the limiting harmonic map 
Q(°): 

Proposition 5 Let C K."^ 6e a simply- connected hounded open set with smooth boundary. Let Q^^''^ G 
W^''^{fl, So) be a sequence of global minimizers for the energy Flq[Q] in the admissible space Aq. Assume 
that as Lk ^ we have Q^^"^ Q(°) in W^-^{n, So). 

Let K d Q be a compact set which contains no singularity of Q^^\ Then 

lim Q^^^^x) = Q^^^x), uniformly for x e K (76) 

k — *oo 



Proof. From the hypothesis and Proposition |3]we have that /s(Q^^'=-'(a;)) uniformly in K. Thus for 
any Eq > there exists a Zq > such that for Lk < Lq we have that (a;) — s+ {n{x) (X) n{x) — ^Id) \ < Eq 

for all X G K (and for each x £ K, we have n{x) G S^). Thus we can apply Lemmas [5j [6]and [T] 

In order to show the uniform convergence it suffices to show that we have uniform (independent of Lk) 
Lipschitz bounds on Q^^''\x) for x G K. We reason similarily to the proof in Proposition H] (see also (Sj). 
We first claim that there exists an ei > so that 



Ve G (0, £i), there exists ro{e) depending only on e, fl, K, and boundary data Qfc so that 



1 



i|VQ(-^'')(a;)p + M^^Mldx < e,Va; G K, provided that Lk < L^{e,ro{e)) (77) 
'"0 JKnBrgix) 2 Lk 
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In order to prove the claim let us first recall that (5*^°^ has no singularities on the compact set K. Thus 
there exists a larger compact set K' with K C K' and a constant C > so that IVQ^^H^^)! ^ C,Vx € K' . 
We choose £i > so that B{x, ei) n K C K' hence for an arbitrary e G (0, ei) there exists ro(£) > so that 

- / i|VQ(°)(:.)pdx<^ 
''0 JKns^(,(K) ^ 

provided that x G K and ro(e) is chosen small enough. We also have, from the W^''^{n, Sq) convergence of 
Q^^*") to (5^°^ that there exists Z(e) so that 

-/ i|VQ(^'=)(x)p(x)dx< - / i|VQ(«)(x)|2da; + |,VLfc<L(e) 

'^0 JKnBrgix) ^ ^0 J/-s-nSrQ(x) ^ 

Recall from the proof of Lemma [3] that limL^^Q J^^ ^'^^^ ^^^^ '^^ ^ ^- Hence there exists Z(e) so that 

-fn ^^^'^Lfc'"'^'^^^ dx < |, VL < L{e). Letting L*(e, ro(e)) — min{Z, L} and combining the two relations above 
we obtain the claim ((77|) . 

Choosing e > smaller than the constant Ci from Lemma [71 we apply Lemma [7] to conclude that 
|yQ(ifc) (^2,^1 Qa^]2 be bounded, independently of Lk, on the set K. The uniform convergence result now follows. 
□ 



4.2 The analysis near the boundary 

In this section we consider the behaviour of a global minimizer Q*-^-* near the boundary, dQ, in the limit 
L ^ 0. For G 911 we define the region fir to be: 

nr''= Qr\Brix^),r>0. (78) 

Lemma 8 Let be a simply- connected, bounded open set with Lipschitz boundary. There exists a constant 
D > 0, depending only on f2, and a constant j'q > such that for all r < r^ and for any x'^ £ dft, we have: 

n^{dnr\Br{x°)) < Dr^. (79) 



Proof. Since fl has Lipschitz boundary, we have that for any G dH., there exists a A(x°) > and an 
orthonormal coordinate system X = (xi, X2, X3) such that a;'' = (0, 0, 0) and there exists a Lipschitz function, 
/j.0 : ^ M, with the property 

U^o {x G n, \xi\ < X{x°),i = 1,2,3} {a- G M^ 2:3 < f^o{xi,X2), \xi\ < A(a;°),i 1,2,3}. 

As il is bounded, it is necessarily uniformly Lipschitz (sec for instance [TO])- Hence, for each a;" G dil, we 
can choose the system of coordinates as before such that there exists a constant c > 0, independent of x", so 
that ||V/^o|| < c,VxO G dn. 

Letting rQ — A we have: 



n Brix")) < / y/1 + \Vf^o{xi,X2)\^dxi dx2 < f \ ^ C? dxi dx2 = 4^1 + 



Vr < Tq. □ 

We have a boundary analogue of the interior mononicity lemma, Lemma [2J namely 



20 



Lemma 9 (boundary monotonicity) Let Q be a simply- connected bounded open set with smooth boundary. 
Let Q'^^ be a global minimizer of Flg[Q] *^ the admissible class Aq. Let 



dV 



Then there exists tq > so that 

d_ 

dr 

where the positive constant C is independent of L. 



Er > {a^,b'',c', Qb, ro,n) , VO < r < ro 



(80) 



(81) 



Proof. Step 1 We assume that the domain fl is star-shaped. Then the proof of (|5T|) closely follows the 
arguments in |20| combined with an idea from [3]. 

Recall that Q'^^' satisfies the equation: 



L 



dQ 



(82) 



In what follows, we drop the superscript L for convenience. 

We multiply both sides of ((5^ by (xp — Xp)Qij^p and integrate over 17,,. Then 



kQij,p{xp - Xp)nk da - jVQp + Qij,kQij,kpixp - x°) dx (83) 



where n is the unit outward normal to dflr and da is the area element on dflr- 

The integral Jq^ Qij,kQij,p{xp — x'^)nk da is evaluated by considering the contributions from dQ,nBr{x'^) 
and r2 n dBr{x^) separately. On n dBr, n{x) = ^^ZxO\ so that 



nndBr 



Qij,kQij,p{-^p Xp^rij^ da 



nndBr- 



dQ 



dr 



da. 



Similarly 



,kQij,p{-^p •^p)^k da 



{x — x'^) ■ n 



dQ 



dn 



+ {X -X^)-T 



dQb dQ 
dr dn 



da 



where t{x) S is the tangential direction to the boundary at x G dil. 
In order to estimate J^j Qij,kQij,kp{xp — x^) dx we note that 



Qij.kQij^kpi-^p -^p) 



d 

dXn 



{xp-xl)-\yQ\^ 



:|VQ|^ 



and therefore 



|VQp + 0,j.fcQ.y,fcp(a-p - x°) dx^ [ {xp ~ xl)\\VQ\^np da - j \ |VQ|' dx. 

The surface integral over d^r can again be expressed in terms of separate contributions from 917 n Br{x^) 
a,ndnndBr{x°). 



Combining the above, we have 

}ij,kk{-^p Xp)QijyP dx 



dx 



il 


dQ 


2 


ivgp 


\ J nndBr- 


dn 




2 



da 



onnBr 



{x — x°) • n 



1 


dQ 


' 1 


dQb 


2' 


2 


dn 


2 


dr 





da 



(84) 
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In we use the fact that |V(5p = 





2 




dQ 


+ 




dn 




dr 



on dft. 



Using the same sort of arguments as above, we compute 



djjB 



d 



{xp - Xp)Qij^p dx ^ — — fB{Q)ixp - Xp) - 3/b((5) dx 



where fsiQ) = fsiQb) = on dfl (from our choice of the boundary condition Qb in (fTTjl ). 
Equating fM]) and ([55)1 we obtain 



(85) 



|vg|^ 



da; + r 



(/ 


dQ 


2 




\JnndBr 


dn 




2 



C?<T 



(x — • n 



1 


9Q 


2 


1 




2' 


dcr + 




dQ 


2 


(9n 




2 


9r 




dn 



nndBr 



fsiQ) , „ 

r acr — 3 

L 



fsiQ) 



dx. 



(86) 



We muhiply both sides of ((86|) by and after some re-arrangement, obtain 



|vg|2 , /b(q) 



+ - 

^ ^ JnndB, 



|VQ|2 , /b(Q) 



(icT = 



1 



JnndB^ 



1 



ag 



9?! 



da 



Ib{Q) 



dx 
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— / (x - xO) • n 



ag 



i9n 



sgt 



dr 



c?<T + — 



, On dQb dQ 
[x - x") ■ r— — • — da. (87) 
OT on 



For a star-shaped domain {x — x^) • n > on dfl. Therefore, the negative contributions to the right hand 
side of ([H7]) are Iou^bM ~ ' 



da and potentially ^ Igs^^g^ix - x") • • || dcr. The first 

integral can be easily estimated since Qb is known. Using the fact that < Cs^ for some C > (as 

Qb G C°°(3il) by hypothesis) where 5+ is defined in (jlOp . we have that 



1 
2^ 



(x-x°) -nl 



dQb 



da < Crs%. 



Here we have used | (x — x") • n| < r and Lemma [5] 
Using Cauchy-Schwarz, we have 



1 



^ JdnnB,. 



dQt 



dr 



1/2 



(88) 



dnnBr- 



dQ 



dn 



1/2 



(89) 



The first integral on the right hand side is easily dealt with i.e. 



annB,. 



dQb 



da < Cs'^r'^, from Lemma [51 



The second integral involving 



is estimated using Lemma [TUl 





dQ 


/ 


dn 



da <G{Qb,n) 



where G > is a constant independent of L. 
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Combining the above we have that 



2 L 



da>-Crsl-G{a^,b^,c'^,n) (90) 



where C and G are positive constants independent of L. We note that 



dr 



L 



^ JnndBr 



|VQP fB{Q) 

2 L 



da 



(91) 



and the above holds for any < r < tq where ro is the constant from Lemma |H1 Therefore 



—£r > ~G' {a^,b^,c'',Qb,ro,n) 

where G > is independent of L. 
Step 2: General domain Q. 

We do not assume that the domain fl is star-shaped and take into account the perturbation terms induced 
by omitting this assumption. As in |20| . the boundary regularity of the domain implies that 



{x — x^) ■ n> \{x — x^) ■ n\ — cr^ 



where c > is independent of r or a;*^ e dfl. Then 



1 
2^ 



{x — x^) ■ n 



dQ 



dr 



da > 



2r2 



I {x — a;°) • n\ 



dQ 


2 p 


dQ 


dn 


^ JdnnBr- 


dn 



da. 



(92) 



(93) 



The inequality ([5T|) now follows from Lemma [TUlD 



Lemma 10 Let Q^^^ be a minimizer of Flq[Q] in Aq (see U2\) ) for a fixed L > 0. Then 

2 

da < G{Qb,VL) 

where G > only depends on the boundary condition Qf, and f2. 





dQ^^^ 


/ 

Jdn 


dn 



(94) 



Proof. The proof follows closely the arguments of Proposition 3 in [3j . Let 1/ : ^ R-^ be a smooth 
vector field on Vl such that y = n on dQ,. We drop the superscript L for convenience. We multiply (|14p by 
VpQij,p and note that 



Qij^kkQij.p^p dx 



1 


dQ 


da~ f 


Jdn 


dn 


Jn 



(95) 



Proceeding similarly as in [3] , we have that 



d 



j,k-Q^ iQij,pVp) dx ~ I Qij,kQij,kpVp + I Qij,kQij,k-^^ dx 



dVr, 



" dXn 



|VQ|^ 



an 



da + (s^) , as L ^ 



Thus, 



Qij^kkQij^p^p dx 



an 



dQ 



an 



da + O (s+) , as L ^ 



(96) 



(97) 
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On the other hand, 

1 
1 



^1 Q V ^I^d. 



T / fB{Q)y -V dx <0{s\), asL^O 



(98) 



since j- dV < C{il)s'^ from energy minimality and fsiQb) = by our choice of Qi, 



Equating ^ and we obtain 



and 



'an 

now follows. □ 



dQ 



an 2 2 





dQ 


2 


dQb 


/ 

Jon 


dn 




dr 



da < C{n)sl 



(99) 



We now prove the uniform convergence of the bulk energy density, /s((5'^^), to its minimal value, on 
compact subsets, K G Q, that do not contain defects of the limiting harmonic map Q^^\ This extends the 
result in Proposition [3] where the uniform convergence is proven only for K C fi. 



Proposition 6 Let Q^^' denote a global minimizer of F]^g[Q] the admissible space Aq defined in (Tl 
Consider a sequence {Q^^''^}ken which converges to a limiting harmonic map Q^^^ strongly in W^'^{n, So) as 
Lk ^ 0. 

Let e dft be a boundary point. We assume that the region f2,. in 1^78}^ contains no singularity of the 
limiting harmonic map Q^^K Then 



Mm fBiQ^^''Hx)) = WxeQr 



(100) 



and the limit is uniform on fJ,. 



Proof. We set a = /s (Q^^'H^")) > 0. Consider the region flp C where p < r < rQ (here is the 
constant from Lemmas |8] and [9|) . Then the boundary monotonicity inequality (jSip implies that 



£p<£r + Cia^,b^,c'',Qb,,n) (r-p) 



(101) 



for p < r < ro- 

Take an arbitrary e > 0. Recall that Q^-^*"^ Q^^^ in W^'"^ as and Qr contains no singularities of 

Q^^\ Using the same arguments as in Proposition |31 we have that there exists an ri < minjro, e} and L > 
( both depending on e) so that for Lk < L 



- I -IVQ^^*-)!' dx<e 



Similarly, we have that there exists an i > (depending on e) so that 



(Lk)] 



dV <e 



for Lk < L { see the proof of Lemma Combining the above, we obtain 

fB (g(^'=' 



1 

P Jn„ 



Lk 



dx <£p< C'e 



(102) 



for any p < ri and for Lk < min{L, L} where the constant C > is independent of Lk ■ 

Using arguments very close to those in [3] ( Lemma A.2 and the way it is used in Step B.l of the proof of 
Theorem 1) together with Proposition [31 one can easily obtain: 



|VQ 



\L=°[n) 



< 



H{a^,b^,c^,n) 



(103) 
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On the other hand, fsiQ) is a Lipschitz function of the Q-tensor and one can infer the foUowing from (|103p 
and Proposition [3l 

ll^^-K^O'l-^^^^ilf^ (104) 

so that 
We take 

P = 



2D{a^,b^,c^,n)' 
There exists a constant so that 

\^p\ > imp' 

(see also [3] for the 2D version of the above) Combining the above with (|105p , we obtain the foUowing 
inequahty 

rfx>— a = p] =1 ^,,, ^, , (106) 



where the constant D > is independent of Lk- Combining p02p and (|106p . wc have that 

< d" {a^,b'^,c^,Qb,n)s (107) 

where D > is independent of Lk- The upper bound p07p is independent of x'^ and e > was chosen 
arbitrarily. Therefore, Proposition [6] now follows. □ 



5 Consequences of the convergence results 



In this section, we discuss some consequences of the convergence results in Propositions HI \E\ and \E\ We 
consider a sequence of global minimizers {Q^^''''}keN converging to a limiting harmonic map Q<-°\ From 
Proposition [SI we have that for a ball B{x, tq) C ft, where B{x, tq) does not contain any singularities of Q'-"' 



)(^'=)(2/)-QW(y) <e(Lfc) yeB{x,ro) 



(108) 



where e(ifc) 0^ as Lk 0. Further, the small energy regularity in Lemma [7] implies that for Lk sufficiently 
small. 



eL, 



Lk 



VQ^^'-Hy) <C(a^6^c^^!) yeB{x,ro) (109) 
where C > is a positive constant independent of Lk- Therefore, for sufficiently small Lk, one has 

|VO(^^)(y)P<C(a^6^c^^7) 

fB{Q^'"'\y))<C{a^b^c^n)Lk yeB{x,ro). (110) 
One immediate consequence of the uniform convergence in (jlOSp and the bounds in piOp is the following 

Lemma 11 Let Q*-^^ denote a global minimizer of Flg[Q] in the admissible class Aq. Consider a sequence 
{Q^^^^}keN which converges to a limiting harmonic map Q^^^ strongly in W'^'^{Q, So) as Lk — > 0. Let x £ fl 
be such that i?(x,ro) H (for tq smaller than the tq used in Lemma 0), does not contain any singularities of 
the limiting map Q(o). Then 

where IS^"^"^ - s+\ < e,{Lk), \R^'"''>\ < t2{Lk) (HI) 
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with n^^''^ 771^^*=-' and p^^''^ unit eigenvectors of Q^^''\ and ei(Lfc), e2(ifc) — > 0^ as Lk — *■ 0. Secondly, if 
X € fl is an interior point such that B{x,rQ) C ft does not contain any singularities of Q^^\ then we also have 
that 

(n(^'=).n("))2>l-e3(ifc) (112) 
where e-i{Lk) as ^ and n^^^ has been defined in [Tl 



Proof The representation piip is a direct consequence of Propositions [21 [4] and [H In the following we 
drop the superscripts Lk for convenience, but keep the superscript in Q^^^ and n'^^\ From Proposition [¥] 
and Proposition [HI we have that 

fB{Q{y))^0 asLk^O 

for y G B{x,rQ) D Q where B{x,ro) D does not contain any singularities of Q'^^\ The bulk energy density 
fsiQ) is a smooth function of the order parameters (S*, R) in Proposition [21 Therefore, as fB{Q{v)) ^ 0, the 
corresponding order parameters {S,R) approach the bulk energy minimum defined by {S,R) — (s-|-,0) and 
the inequalities piip follow. Further, if B{x,rQ) C fl, then the uniform convergence ()108p holds. A direct 
computation shows that for 



Q = 5 ( n ® n — -Id ) + R{m ® m — p ® p) , 



we have 



The lower bound on 



now follows from (jlOSp and the fact that \S — s+| < ei(Lfe), \R\ < e2{Lk) for 



sufficiently small values of i^. □ 

Proposition 7 Let Q^^^ denote a global minimizer of Flg[Q] the admissible space Aq. Consider a se- 
quence {Q'^''^}a;gn which converges to a limiting harmonic map g(0) strongly in W^^'^{n,SQ) as Lk 0. 
Then Q^^''^ converges uniformly to the limiting harmonic map Q^^\ away from the singular set of Q^^\ in 
the interior ofQ. Let K <Z Vl be an interior subset that does not contain any singularities of Q^^\ Then (i) 



/3(Q(^'')(y))<CLfc 



<D^Lk yeK 



(113) 



where C and D are positive constants independent of Lk. (ii)(rate of convergence of eigenvalues) Let |Aj^'°^| 
denote the set of eigenvalues of Q^^''^ and {Xi} denote the set of eigenvalues of Q^^\ Then 



Af''^(tj)-A.(y) <a(a^5^c^)Lfe yeK;i = 1...3 



(114) 



where a is a positive constant independent of Lk ■ 



Proof (i) This follows directly from piOp and Proposition [H In Proposition we obtain a lower bound for 
fsiQ) in terms of |(5| and f3 and in piOp we have an upper bound for fsiQ) in terms of Lk as shown below 



101- 



6V6 



m)\Q\' < Ib {Q{y)) < C {a\b\c\ n)Lk yeK 



(ii) In the following we drop the superscripts "(Lfc)" for convenience. From piop . we have the following 
upper bound for the bulk energy density on the set if C O 



fB{Qiy))<CLk yeK 
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where C is a positive constant independent of Lk- Using the representation formula we have that 

Q = S [n(E)n-^Id] + R{m(g)m~p(E)p) (115) 



where 

\S - s+\ < e^iLk) = 0(1) (116) 

and 

<e5(ifc) =o(l). (117) 

A direct computation shows that 

2 2S^ 
QijQij = -^S +2R QipQpjQij ~ 2SR . (118) 

From (|116p and (|117p . we represent Q on the subset X C as fohows: 

Q = (s+ + e) (K) n — ^Id^ + ^{m®m— p®p) (119) 
where |e|, I7I = o(l). Using pisp . wc find that 

|Q|' = ^(4+e' + 2s+e)+272 



and from the maximum principle (Proposition [3|) , 



\Q{x)\^<\sl xeK. 

This necessarily implies that e < 0. 

The bulk energy density fs is given by 

fBiQ) = y (4 - S') + ^ {si -S')-j (4 - 5') - + ^SR^ + ^-^S'R' + c^R' (120) 

where we have merely expressed trQ^ and trQ'^ in terms of the order parameters S and R. We write the bulk 
energy density as the sum of two contributions - 

fB{Q) = F{S) + G{S,R) (121) 

where 

Fis) = ^ (4 - + ^ (4 -s')-j (4 - s') 

and 

052 2c^ 

G [S, R) = -a^R^ + —SR^ + —S^R^ + c^R'^. 

o o 

The function F{S) is analyzed in (|149p : the function F{S) is bounded from below by 

F(5) >L»(a^6^c2) (S'-s+)^ D{a^,b^,c^)>0 (122) 

Similarly, since 2c?s'\_ = 6^s+ + 3a^ and < s+ — 5* = o(l) (for sufficiently small), we have the following 
inequality 

o o 252 , 2c2 , , 52 

^a^R^ + S'i?2 + S^R^> (123) 

3 3 2 



Combining ([T^ . ([Hg]) and (fTTUll . we obtain the following 



i?(a2, c^) £2 + ^s+j' + 0^74 < {Q{x)) <C Lk (124) 
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from which we deduce 

< Ci Lk and 7^ < C2 Lu, 
where Ci,C2 are positive constants independent of Lk- The inequaUties ()114p now follow. □ 
Next, we have a lemma about the leading eigenvector n in the representation (jllip . 



Lemma 12 Let Q = S {n® n — \ld) + R (m ® m — p ® p) with S > 8\R\ and n,m,p £ pairwise perpen- 
dicular. Then the minimum of 



a(E) a — -Id 



with a S is attained by a = ±n. 



Proof. A direct computation shows that 

2 r. 



Q - s+ ( a (g) a — i/d 



2S'.s+ (n • a) — 2s^R ((m-a) — (p ■ a) 



2R^ 



2s+R - 2s+ (5 + R){n- ay - As+R (m • ay 



(125) 



where in the last line of p25p . we use the equality (n • a)^ + (m • a)^ + {p ■ 



1. Since S > 8\R\, one can 



immediately verify that (jl25p is minimized for (n • a)^ = 1 or equivalcntly a = ±n. □ 

We can now provide a result about the regularity the leading "eigendircction" n^n G M^^^ where n S §^ 
is the leading eigenvector. For a thorough discussion about the relationships between the regularity of the 
eigenvector n e §^ and that of the eigendircction n® n £ M^'^^ see [5]. 



Corollary 1 Let Q'^' denote a global minimizer of Fijg[Q\ ^n- the admissible class Aq. Consider a sequence 
|Q(ife)|^gj^ lufiicfi converges to a limiting harmonic map Q^^^ strongly in W^'^{n, Sq) as Lk 0. Let K <Z Vl be 
a compact subset of fl that does not contain singularities of the limiting map Q^^\ Then, for L^ small enough 
(depending on K), Q^^''^ can be represented as in Ull\) on the set K G ft and the leading eigendircction 



Proof From piip . we can represent Q'^*"' as 

where |S'(^'') - s+| = o(l), o(l), and n^^"', m(^'=),p(^'=) e ^rc the eigenvectors of Q'-^'^K 

Let 7r(Q) be the nearest neighbor projection onto the manifold of global minimizers of the bulk energy 
density, denoted by Qmin = {s+ (a® a — , a S S^} as in Namely, n{Q) associates with each Q' , (in 
a neighborhood of the manifold Qmin) an element Q* G Qmin such that 

|Q'-Q*|= min \Q' - Q\ . 

Q£ Qmin 

The projection tt is defined only in a neighborhood of the manifold Qmin and moreover tt{Q') G C°° {Sq, Qmin) 
(see, for instance, [S]). The Lemma fTDand Lemma fT^ show that in our case 



7r(Q 



= s+ [ n'-^"^ (g) n^^"'^ - -Id 



Therefore, the tensor 



{n'^^''^ (E) n'^^'' 



-Id) G {K, So) , 



(since s+ is a constant) and the conclusion of the lemma now follows. □ 
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6 Biaxiality and uniaxiality 
6.1 The bulk energy density 

Our first proposition concerns the stationary points of the bulk energy density. 
Proposition 8 \21}j Consider the bulk energy density fsiQ) given by 

fB{Q) = -\trQ' - \trq^ + \ {trq-f + ^4 + ^4 - ^4. (126) 



Then fsiQ) attains its minimum for uniaxial Q-tensors of the form 

1 
3 

where 



s+ (n(g)n - I ] , (127) 



62 + VFT24^ , 
= ^ ^^^^) 

and n : ri — > S*^ is a unit eigenvector of Q . 

Proof. Proposition[8]has been proven in [21] and we reproduce the proof in the Appendix for completeness. 

□ 

In the following proposition, we estimate fsiQ) in terms of |(5| and the biaxiality parameter (3(Q). 
Proposition 9 Let Q Sq. Then the bulk energy density fsiQ) is bounded from below by 

MQ) > 



/ \ 2 



|Q|-\/|s+j +-^/3(0)|g|3 (129) 



where S-|_ has been defined in U28\) . 

Proof. From Lemma [l] we have the inequality, 



From the definition of fsiQ) and s+ in (|126p and (1128p . we can obtain a lower bound for fsiQ) in terms 
of IQI and f3{Q) as follows i.e. 



pmQ\'- (131) 



> 



6^/6 



The bracketed term in p3ip can be further simplified by carrying out a series of calculations. Consider 
the function 

r)2 ?,2 ^2 

^ ^ 2 3V6 4 ^ ' 

The stationary points of f{u) are solutions of the algebraic equation 

V6 



/■ (u) = It ( c^it^ - -a^] =0 (133) 
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and one can readily verify that f{u) attains its minimum for 

[2 



Umin = \l gS+ 



(134) 



The bracketed term in (|13ip is non-negative by virtue of (|132p - ()134p . Further, let S = \Q\ — where 



„22 2 _ b_ 

^ 3 + ^ Vey 3 



s_|_ + by the definition of s+. Then 



= 5 



2 2^2 2 s 
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.2 2 

s+ + s.c 

2 3 + + 



(135) 



(136) 



The coefficient of 5 vanishes from the definition of s-|_ in (jl28p . The coefficients of 5 and 5 are positive since 



> 



262 



3^ ^"^ 3^/6 ' 3^/6 



We substitute p77)) into to obtain 



3V6 
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> 



|c2s^ + a2 



(137) 



(138) 



and on combining (jl38p witli (|13ip . tlic lower bound (|129p follows. □ 

The bulk energy density, /s(Q), can be equivalently expressed in terms of the order parameters s and 
in Proposition [11 as shown below 

Proposition 10 Let Q ^ Sq be represented as in Proposition]^ 



s [ n (E) n — 2^*^ j + ^ f™ ^ — -^Id 



with either < r < ^ or|<r<0. Case (i) Non-negative order parameters, < r < ^ with < s < s+; 
where S-|_ is defined in il28\) . Then the bulk energy density, fsiQ), is bounded from below by 



fsiQ) > (.s+-s)S(a^6^c2) + 



r(s — r) 



9 ^ ' 27 

where 7 {a^ ,b^ ,c^^ is an explicitly computable positive constant. 

Case (a) Non-negative order parameters, < r < f and s > Then 



562 

'Sa^ + b^s~2c^s^) + —r^s < s < s+ (139) 



fsiQ) > 



min {^{s- s+f , ^ (VSs - 2s+ 



^52 „3 ir^s-rf 



(140) 



where r is an explicitly computable positive constant, independent of a ,b ,c . 
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Case (Hi) If ^ < r <0, then 

fB{Q) = fsi-Q) + ^ (2|sp + 2|r|3 - isM - Mr^) , (141) 

where —Q G Sq has positive order parameters < — ?' < — f fsi^Q) can be estimated using il39\) and 
mO^ - In particular, 

.3 ^2 



MQ)>-i^-f ^9-y^+j>o (142) 

for Q-tensors with ^ < r < Q. 

Proof. From Proposition [U it suffices to consider the two cases < r < | and | < r < 0. 

Case (i): We can explicitly express the bulk energy density. fsiQ), in terms of s and r as follows - 

fsiQ) = -y (s' +r^-sr)- — (25^ + - 3s\ - 3sr^) + 



c 



262 



+ — (s^ + r^ + 3s^r^ - 2sr^ - 2s^r) + — s| + —s% - —3%, (143) 



9 ^ ^3+27 
where we have expressed trQ^ and trQ'' in terms of s and r 

2 

trQ^ = - (s^ + - sr) 
3 

and 

trQ^ = - (2s^ + 2r^ - 3sV - 3sr2) 
The function fsiQ) consists of two components - 



fB{Q) = Fis) + Gis,r) where 

"27 ^" ' 9 



0(5, r) = y (sr - r^) + ^ (Ss^r + 3sr2 - 2r^) + — (-2sV + 3s^r^ - 2sr^ + r^) . (144) 

Recalling that 2c'^s'^ = b^s+ + (from the definition of s_)- in (|128p ). the fimction F{s) can be expressed in 
terms of(5 = s+ — s>Oas follows - 

F{s) = — {l8a^s+ + 66^4 - 12c^4) + 
The coefficient of 5 vanishes by virtue of the definition of s+ in (|128p . We note that the function 
attains a minimum for 

(5„,in = 2s+ - — > s+ (147) 



and, therefore. 



G(5) > G(s+) = 1 (262,s+ - 9c24) . (148) 



We substitute (|148p into (|145p to obtain the following lower bound for F{s) 



F(s) > (s+ - s)' . (149) 
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We can analyze the function G(s,r), in ()144p . in an analogous manner. Let 7 = f e [0, i] . Then 



G{s, r) = 



2c' 



s s 

3 27 9 



I 2 2 

+ 7 s 



35^ 
27 



+ 7^s^ 



27 



9 +^^^150) 



The coefficient of 7 is non-negative for all s < s+. Using the inequality 7 < ^, one readily obtains the 
following lower bound for G(s, r) - 



G{s,r) > 75^ 



362 



2c' 



s s 

3 27 9 



2 2 

+ 7 s 



262 



2c2 



> 



> 



562 



±^{3a' + b\s-2c's') + ^-^^'^ 



(151) 



Combining (|f 49p and p5f p . the lower bound for < s < s_|- in p39p follows. 

Case(ii) The case s > s+ can be dealt with similarly. For any Q E Sq with < r < §, we have that 



(152) 



For s > s 



+ 1 



f > ^ and 



/5(Q) > V 



r'{s — r)" 



(153) 



where P{Q) is the biaxiality parameter defined in ([22|) and ?/ is a positive constant independent of a2, 62 or c2 
or L. Combining (|f 52p . (|f 53p and (|129p . we readily obtain the lower bound 



fsiQ) > 



2 2 2 I 2 



> 



2 



\Q\ 
2 



6^/6 



/3(Q)IQI'> 



mm 2 (s - s+f , - iVSs ~ 2s 



7-2 (s — r)^ 



(154) 



where r is an explicitly computable positive constant. 

Case (iii) Finally, we consider Q E Sq with negative order parameters § < r < 0. In this case, one can 
directly check that 

trg^ = 7: (2s-'' + 2r^ - 3s2r - 3sr2) < 



9 



and therefore. 



« -12 ^ ,..^3 



/i3(Q) = -ylOI'-ytrg^ + jlQI' 



a2 2 262 ^ ^2 ^ 



a 

y 



,9 ^ . c2,^,, a2 „ 262 2 262 I 

P - ytr {^Qf + -|Q|4 + -4 + —4 _ _4 + _ jtrQ^ 



(155) 



since ^tr {-Qf = --ytrQ^ and --ytrQ^ = ^ |trQ^|. The inequality (|14ip follows from (|155p upon express- 
ing trQ"^ in terms of s and r. 

For p42p . it suffices to note that for s,r < 0, tiQ^ < and therefore. 



fsiQ) > -ylQI' + jlWr + ys:; 



262 
'27 



.3 /^2 ^2 



r — sr Sa 

' 3 V 9 3 ^ 
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A straightforward computation shows that the function 

2 2 



4 

2 a 



and 

c3 ^2 



The inequality (|142p now follows. □ 



3 V 9 3"^+j ^ 4c2' 



Remark 2 One can readily obtain lower bounds for fsiQ) in terms of the order parameters (S, R) in Propo- 
sition [H following the methods outlined in Proposition \1U[ The details are omitted here for brevity. 

Remark 3 Relation shows that if fB{Q^'"''\x)) as ^ then Q^^'-'' cannot have an (s,r) 

representation with | < r < 0, if Lk is sufficiently small. 

In view of Propositions [3] and [HI we can make qualitative predictions about the size of regions where a 
global Landau-De Gennes minimizer Q* can have \Q*\ << ^/|s+ and the size of regions where Q* can be 
strongly biaxial. 

Proposition 11 Let Q* be a global minimizer of Flg[Q] in I113\} . in the admissible class Aq defined in il2\) . 
Let 17* = |a; e fi; \Q*{x)\ < ^^J^s+Y Then 

\n*\<a^-^—-f\Vn^°\x)\^dx, (156) 



[c^sl + a^) Jn 

where n'"-' is defined in hl6\) and a is an explicitly computable positive constant independent of , 6^, or L. 

Proof. From Proposition wc have that 

/b(Q*(^))> -(c'4+a')4, xen* (157) 
a 



for some explicitly computable positive constant a, since \Q* \ < iy^|s+ = on Q* . On the other hand, 

recalling the definition of Q'"'' in (fT5|) and since Q* is a global minimizer of Flg[Q], wc have that 

/ /i3(g*(x)) dx < TlgIQ^'^^ = I .hiQ^''^) + ^|VQ("f dx = Lsl [ |Vn(")|2 dx, (158) 

since fsiQ'^^^) = everywhere in Vl. Substituting (|157p into p58p . wc obtain 

- {c^sl + a") sl\n*\ < Lsl I dx, (159) 

from which the inequality p56p follows. □ 



Proposition 12 Let Q* be a global minimizer of Flq[Q] in (f in the admissible class Aq defined in 112\) . 
Let = |.T G ri; \Q* {x)\ > ^<fj^s+, P{Q{x)) > a| for some positive constant X. Then, 

\^^\<a-^ I |Vn("'|2dx (160) 
'^s+o Jo. 

where n'*'' is defined in Mb]) and a is an explicitly computable positive constant independent of a? , b^, or L. 
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Proof. From Proposition [9l we have that 



fB{Q*{x)) > l^/3{Q*{x))\Q*{x)\^ > -b'Xsl xen^ (161) 
6V6 a 



for some explicitly computable positive constant a, since \Q*\ > |s+ — '^^+ the other hand, 

recalling the definition of Q^^\ (fT5|l . and since Q* is a global minimizcr of FlqIQ], wc have that 

/ fBiQ*ix)) dx< [ /B(0^")) + ^|Vg(°f dx = L4 / iVnWpdx (162) 
since fsiQ''^^) = everywhere in r2. Substituting (|161|1 into p62|) . we obtain 

-b^As^ll^^l < i4 / iVn^p dx, (163) 

from which the inequality (|160p follows. □ 

Proposition 1 1 1 1 is relevant to the size of defect cores in global energy minimizers whereas Proposition 1 1 21 is 
relevant to the equilibrium behaviour far away from the defect cores. 

6.2 Analyticity and uniaxiality 

We define a new biaxiality parameter (3{Q) as follows: 

/3(g) =^tr(g2))3_6(tr(Q3))2. 

Then /3(g) > with /3(g) = if and only if Q is uniaxial i.e. g = s (n (g) n - ^Id) for some s G M \ {0}, n e §^ 
or g = 0. The function /3(g) is a real analytic function of Q and this is particularly important given that 
global energy minimizers of the functional (subject to smooth boundary conditions) are real analytic: 

Proposition 13 Let il be a simply-connected bounded open set. Let g'^-* be a global energy minimizer of 
Flg[Q] in hl3\) in the admissible space Aq. Then g^^-* is real analytic in fi. 

Proof. 

We drop the superscript L from Q^^^ for convenience. As — ^tr {Q^) — ^-tr {Q'^) + % (trg^)" is bounded 
from below (sec also the Appendix) we have that there exists an global energy minimizer satisfying the 
Euler-Lagrangc system: 

LAg,y = -a^gy - b" {Q,kQk3 - y tr(g2) j + c^g,, tr(g2) 

For g an solution of the equation one uses ^ L^ (in M'^) and Holder's inequality to obtain that 
the right hand side of each equation is in L^. Elliptic regularity gives that Q G ^ W^'^ ^ L°° hence the 
right hand side of the equation is in . Elliptic regularity gives Q e and one can continue bootstrapping 
to obtain the full regularity allowed by the regularity of boundary data and that of the domain. 

In order to prove the analyticity we use a general abstract result due to A.Fricdman.[12j. We define growth 
classes as follows: let Mn be a sequence of positive numbers. Then a function F : C°°{D) — *■ C, with D C M"^ 
an open set, belongs to the class C{Mn; D} if for any closed subset C D there exist constants FIq, H with 

\d^F{x)\ < HoWMj,x e Do 

where we have used multiindex-notation [d^ F ~ d-'i . . .d^/F;j = 'Ef^iji). Let us observe that C{n\;D} is 
the class of functions analytic in D. 

In [12| the following theorem is proved for general elliptic systems: 
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Proposition 14 ( U2}/ .p.4-5) Let u{x) be a real solution of the elliptic system 

$/(a;;u,Vii,V2u, ...,V2'"u) ^ 0,x e Q (ZR'^;u e , I ^ 1, . . . , N 

in il <Z R''. Let E be some open set containing Ei {ii(a;), Vw(x), . . . , V^™w(a;);a; € f2}. Assume that: 
(i) (z C{Mn; X E} and that the A/„ satisfy the monotonicity conditions 
(a) {")MiM„-i < AMn] < i < n, n e N for some A>0. 

Ifue C2'"+"(r2), < a < 1 then u e C{M„_2™+i; n} (where AI-, ^IforieN) 

In our case, for the system ([T4| we have m = 1 and is analytic hence of class C{n!; fl}. The constants 
Mn = n! satisfy the monotonicity conditions (ii) in the theorem, with A = 1. We have that Q G C°°{n) and 
hence by the theorem Q is in the class C{(n — 1)!; SI} therefore real analytic. □ 

Proposition 15 Let Q be a real analytic function Q : fl C M.^ ^ Sq. Then the set where Q is uniaxial or 
isotropic is either the whole of fl or has zero Lebesgue measure. 

Proof. If there is no a; € S7 such that (3{Q{x)) ^ then Q is uniaxial or isotropic everywhere. If there 
exists a P £ n such that j3{Q{P)) ^ then let us consider the lines passing through P. The restriction of Q 
to any such line is real analytic and then so is (3{Q). Thus (3{Q) has at most countably many zeroes on such 
a line. We claim that this implies that the set of zeroes of (3{Q) in SI is of measure zero. 

We assume, without loss of generality, that P = 0. We denote N* ^ N \ {0} and decompose R" n O = 

U„GN- (^Bi \ B 1^^" nnju (^U„gN* (^Bn+i \ Bn r^ Vlj^ . We claim that for any n, i e N* the set {j3{Q)^ (0)n 

{vt n Bi_ \ B_L^ is a set of measure zero. This implies that P{Q)~^{Q) n f2, which is a countable union of 
sets as before, is also a set of measure zero. 

We consider the bi-Lipschitz functions 
1 1, 



/„ : [ -, -] X [0,7r] X • •• X [0,7r] x[0,27r) Bi \ B^,yn, - £ N 

n + 1 n ■• ^ ' " "+1 n 

times 

that realize the change of coordinates from polar to usual cartesian coordinates. 

We have that /"^ (/3(Q)-i(0) n n Bi \ B_i_}j C [;^, x [0, tt] x • • • x [0, tt] x [0, 2tt). We recall that 

n—2 times 

the Lebesgue measure on the n-dimensional product space [;;;qrrj^] x [0,7r] x ••• x [0, tt] x [0, 27r) is the 

n—2 times 

completion of the product measure /^i x ^2 where /xi is the 1 dimensional Lebesgue measure on ^] 

and ^2 is the n — 1 dimensional Lebesgue measure on [0,7r] x • • • x [0, tt] x[0, 27r). Then for any set E C 

" V ' 

n — 2 times 

[;^, i] X JO, tt] X ■ ■ ■ X [0, tt] X [0, 27r) we have 

n — 2 times 



{fxi X fi2)iE) = Mi(i?")M2(rfy) 

'[0, tt] X • • • X [0, tt] x[0,27r) 



where Ev = {x £ i], {x,yi, . . . ,?/n-i) £ E} C [:^, In our case, letting 

E =^ f-' (mr'io) n nnBl\BT:) 
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we have that is made of finitely many points for almost all y G [0, tt] x • • • x [0, tt] x [0, 2tt) (as a consequence 

n—2 times 

of the first paragraph in this proof; because is just the set of the distances to P of the uniaxial or isotropic 
points that are in n Bi_ \ B^_, on a a segment through P, segment that has in polar coordinates the 
direction y G [0, tt] x • • • x [0, tt] x [0, 27r) ). Thus fJ-iiE^) = 0, /i2 — a.e. y hence /ii x H2{E) = thus ^i{E) = 0. 

n—2 times 

As bi-Lipschitz functions carry sets of measure zero into sets of measure zero we have that (3{Q)^^{Q) D 
H. n Bi_ \ is a set of measure zero. On the other hand P{Q)^^{Q) n is a countable union of sets as 

before, hence it has measure zero. □ 

Corollary 2 Let Q(^) be a global minimizer of Flg[Q] *^ the admissible class Aq. Then there exists a set of 
measure zero, possibly empty, VIq in 51 such that the eigenvectors of Q*-^' are smooth at all points x G 51 \ fio- 
The uniaxial-biaxial interfaces, isotropic- uniaxial or isotropic-biaxial interfaces are contained in JIq- 

Proof. The global minimizer 

Q{L) g (jco (f^;_4). The eigenvectors of Q^-^) have the same degree of regularity 
as Q^^^ on sets K C where Q*^'^-* has the same number of distinct eigenvalues i.e. where Q'^^^ is either 
biaxial or uniaxial or isotropic, [26], but not necessarily otherwise |17| . If Q*-^-* is uniaxial everywhere then 
flo = 0. If Q^'^^ is either uniaxial or isotropic on the whole of (i.e. /3(Q^-^^) = in fl), with Q^-^^ ^ 
at some point in il, then let jl = G ^,Q^^\x) = 0} denote the zero-set of Q'^'. Let us observe that 
^ " (IQP) (0) and IQp is an analytic function. By an argument similar to the proof of Proposition [T51 

def ~ 

and since Q{x) ^ for at least one point x G fl, we have that SI has measure zero and we take JIq = 51. 

If Q'^' is biaxial somewhere then Proposition [T51 shows that the set of points where l3{Q) = has measure 
zero. We denote this set by SIq and observe that f2 \ JIq is an open set and the eigenvectors have the same 
regularity as Q^^^ on fl \ flo, see O 
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Appendix 

Proposition 16 \21}/ Consider the bulk energy density fsiQ) given by 



fBiQ)=-Yt^'-jt^''+jM'y- (164) 



Then fsiQ) attains its minimum for uniaxial Q-tensors of the form 



g = s+(n®n-i), (165) 
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where n : fl is a unit eigenvector of Q and 



9 + 

" 4^5 ■ ^ ^ 

Proof. Proposition fTCl has been proven in [21]. We reproduce the proof here for completeness. 
We recall that for a symmetric, traceless matrix Q of the form 

3 



=5: A, 



trQ" = J2^=i subject to the tracelessness condition so that the bulk energy density fs in (|164p only 
depends on the eigenvalues Ai, A2 and A3. Then the stationary points of the bulk energy density Jb are given 
by the stationary points of the function / : M'^ ^ M defined by 

3 3 / 3 \ ^ 3 

/(Ai,A„A3) = -^^A^-^5:A? + ^ 5:An -2<5^A.. (167) 

1=1 i=l \i=l / i=l 

where we have recast Jb in terms of the eigenvalues and introduced a Lagrange multiplier 6 for the tracelessness 
condition. 

The equilibrium equations are given by a system of three algebraic equations 

^■^ 4^ -a^A, -b^A^ ( ^AH A, = 25 for i = 1 . . . 3, (168) 



dX 

or equivalently 



(A, -A,) 



3 

2 

fe 

fe=i 



1 < i < j < 3. (169) 



Let {Ai} be a solution of the system (|168p with three distinct eigenvalues A^ 7^ A2 7^ A3. We consider equation 
()169p for the pairs (Ai, A2) and (Ai, A3). This yields two equations 



(Ai+A2) + c2^A^ = 
fe=i 
3 

-a^ ^b^ (Xi+X3) + c''Y,^l^^ (170) 

k=l 

from which we obtain 

- (A2 - A3) = 0, (171) 

contradicting our initial hypothesis A2 7^ A3. We, thus, conclude that a stationary point of the bulk energy 
density must have at least two equal eigenvalues and therefore correspond to either a uniaxial or isotropic 
liquid crystal state. 

We consider an arbitrary uniaxial state given by (Ai,A2,A3) = (T'~f'~f) ^^^'^ corresponding Q- 
tensor is Q = s (ei ® ei — ■^Id)- The function fs is then a quartic polynomial in the order parameter s 
ie. ^ 

fB{s)^^{-9a^-2bh + 3ch^) (172) 
and the stationary points are solutions of the algebraic equation = 0, 

% = - (-IBa^s - 6bh^ + Uch') = 0. (173) 
ds 27 ^ ^ 
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The cubic equation ()173p admits three solutions; 



s = and s± = (174) 

where 

/b(0) = and < /b(s_) < 0. (175) 

Symmetry considerations show that we obtain the same set of stationary points for the remaining two uniaxial 
choices. The global minimizer is, therefore, a uniaxial Q-tensor of the form 

Q = s+(^®n- i/d^ , n e §2 (176) 

where s+ has been defined in (|166p .D 

Lemma 13 Let Q d Sq. We define the biaxiality parameter P{Q) to be 

/3(g)^l-6 ^^'^'\ . (177) 

(i) The biaxiality parameter (3{Q) £ [0, 1] and p{Q) = if and only if Q is purely uniaxial i.e. if Q is of the 
form, Q = s{n®n— \ld) for some s S IR,n € S^. (ii) The biaxiality parameter, (3{Q), can be bounded in 
terms of the ratio ^, where (s, r) are the scalar order parameters in Proposition]^ . These bounds are given 
by 



ill 



1 - \/]3) < J < i (l + 'Ji^^M . (178) 



Equivalently, 



i+^/l-^ S 3-^1-^ 

where {S, R) are the order parameters in Proposition]^ Further f3{Q) = I if and only if r = ^ or if and only 
if ^ = |. (iii)For an arbitrary Q G Sq, we have that 



Proof: (i) The quantity P{Q) is known as the biaxiality parameter in the liquid crystal literature [23] and it 
is well-known that (3{Q) G [0, 1]. We present a simple proof here for completeness. 

Following Proposition [U we represent an arbitrary Q e 5o as 

-Id] Q <r <- or - <r <Q. (181) 



Since S-j^^^^yg- > 0, the inequality (3{Q) < 1 is trivial. To show [3{Q) > 0, we use the representation p8ip to 
express trQ'^ and trQ^ in terms of the order parameters s and r. 

trQ^ = i (25^ + 2r^ - 3sV - 3sr^) 

trQ2 = ^ (s2 _^ ^2 _ (,g2) 
A straightforward calculation shows that 

{trQ^Y ^ — (45^ + Ar^ - 12sV - 12sr'^ + 26s - 35^^ - 3s^r'^) 
81 
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and 



2 

-=4(1-/3). (184) 



One can then directly verify that 

(trQ2)3 _ 6 (trg3)' = 2s\'' {s - rf > (183) 

as required. It follows immediately from (|183p that /3(Q) = if and only if either s ~ Q,r ~ or s = r. From 
P8ip . the three cases, s = 0, r = and s = r, correspond to uniaxial nematic states (in fact all uniaxial states 
can be described by one of these three conditions) and therefore, f3{Q) = if and only if Q is uniaxial. 

(ii) From Proposition [l] it suffices to consider Q-tensors with either < r < | or | < r < 0. Let 7 = ^, 
then 7 g [0, ^] for the two cases under consideration. The biaxiality parameter, (3{Q), can be expressed in 
terms of the ratio 7 as follows 

(2 - 37 - 37^ + 273) 

(1-7 + 7^)' 
From p83p . we have that 

(2 - 37 - 37^ + 27^)^ = 4 (1 - 7 + 7^)^ - 277^ (1 - jf , (185) 
which in turn, yields the following equality 

'jP^'V. (186, 

(1-7 + 7^) 

Noting that for 7 = [O, ^] , the polynomial 1 — 7 + 7^ > |, we obtain the following upper bound 

(3 < I672 (1 - 7)^ (187) 

and the bounds (|178p readily follow from ()187p . 

One can readily see from (|186p that /3((5) = 1 if and only if ^ = i . The bounds p79p follow directly from 
P78p on noting that 

r = 2R and s ^ S + R. 

One can see directly from (|179p that if /3 = 1, then ■§ = ■!■ On the other hand, if ^ = i, then 7 = ^ and 
P86p implies that (3{Q) = 1. The claims in (ii) now follow. 

(iii) From the definition of the biaxiality parameter in (|177p . we necessarily have that 

trQ3 = ±i^^l_;3(g). (188) 



It is easily checked that 



VW<1-^ (189) 



The bounds pSU]) follow from combining p^H]) and ^TM^ .□ 
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